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Abstract 

We investigate the existence and location of the surface phase known as the "Surface- Attached 
Globule" (SAG) conjectured previously to exist in lattice models of three-dimensional polymers 
when they are attached to a wall that has a short range potential. The bulk phase, where the 
attractive intra-polymer interactions are strong enough to cause a collapse of the polymer into 
a liquid-like globule and the wall either has weak attractive or repulsive interactions, is usually 
denoted Desorbed-Collapsed or DC. Recently this DC phase was conjectured to harbour two 
surface phases separated by a boundary where the bulk free energy is analytic while the surface 
free energy is singular. The surface phase for more attractive values of the wall interaction is 
the SAG phase. We discuss more fully the properties of this proposed surface phase and provide 
Monte Carlo evidence for self-avoiding walks up to length 256 that this surface phase most likely 
does exist. Importantly, we discuss alternatives for the surface phase boundary. In particular, 
we conclude that this boundary may lie along the zero wall interaction line and the bulk phase 
boundaries rather than any new phase boundary curve. 

The phase transitions of a single isolated polymer in solution continue to attract attention as 
single polymers are fundamental components in more complicated modelling scenarios and these 
transitions are not yet fully understood. The collapse transition [1] mediated by the intra-polymer 
attractive interactions and the adsorption transition [2] when a polymer is attached to a sticky wall, 
are two of the key transitions that have been well studied. The situation when both transitions can 
occur in the same system has been studied less intensely due to the difficulty of numerical work 
such as Monte Carlo simulations when two independent parameters compete. However, with the 
advent of more powerful computers and sophisticated algorithms this situation has received some 
attention. 

The standard lattice model for polymers is the self-avoiding walk (SAW) with the sites of the 
walk called monomers. The intra-polymer attraction is modelled by a potential energy Sp associated 
with monomers that are nearest-neighbours on the lattice. These instances will be referred to as 
(nearest-neighbour) contacts. To consider adsorption a wall is introduced, so that the polymer 
is restricted to one side of the wall, or may visit the wall, with one end of the polymer (SAW) 
attached to the wall. Any monomer which visits the wall, other that the one fixed on the wall, is 
given a potential energy £w Tliese monomers will be referred, to as visits. Vcirious different phases 



are conjectured to exist at different values of the ratio of the two energies and the temperature 
(see below). We shall restrict our discussion to three dimensions, as this is where more complex 
behaviour can occur (rather than two dimensions, where more is known due to exactly solved models 
and conformal field theory). 

Consider the simple cubic lattice and the half space z > 0. Consider an n-step self-avoiding 
walk (f with one end fixed at the origin. We define two Boltzmann weights ujp = ^^-^p/^bT 

= e~^'^I^BT ^}igpg is Boltzmann's constant and T is the temperature. Therefore, when 

> 1, = 1 and uJw < 1 the interaction of the walk with the wall is attractive, neutral, and 
repulsive, respectively. When cjp > 1, tjp = 1 and ujp <\ the intra-polymer interaction is attractive, 
neutral, and repulsive, respectively. 

The partition Z„ (01^,0;^) of the model (see Figured]) is defined as 

Z„(a;p,u;^)= ^ J^^^'''^ (1) 

where rj„ is the set of all n-step SAW restricted to the half space z > 0, mp{(p) is the number of con- 
tacts in the walk (p, and m^((/9) is the number of wall visits of the walk 93. The hulk thermodynamic 
(reduced) free energy is given by the limit 

fb{ujp,ujw) = - lim -log(Zn{u}p,u;w)) ■ (2) 




Fi gure 1: A 9-step SAW^ {ji — 9) attached to a wall with three visits (?7i^ 
nearest neighbour contacts (m^ = 2, shown as double helices). 



3) to the wall and two 



The hulk phase diagram is determined by the analytic structure of fhi^^i 



This phase 

structure was studied in a series of papers [3l [H [5] where a schematic phase diagram was proposed 
(see Figure [2]). To discuss this diagram it is worth considering two order parameters. Firstly let us 
define the "internal density" pp as 



Pp = lim Pp{n) = lim 



n 



where Rg (n) is the radius of gyration. Secondly we define the density of wall visits pw as 



lim pyj{n) = lim 



n 



(3) 



(4) 



where (m^) is the ensemble average of the number of visits to the wall. 

For any ratio of wall to intra-polymer energies at high temperatures, the polymer is in the 
excluded volume state and is entropically repulsed from the wall. It is expected that Rg{n) ~ ti^ 
as n ^ 00, where the three-dimensional excluded volume value has been estimated as = f 3 ~ 
0.5874(2) [3. The average number of visits is expected to behave as (m^,) = o(n) (in fact, it is 



expected to be bounded): this has been numerically verified. Hence for uj^ 



1 we have 
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Figure 2: Schematic bulk phase diagram for collapsing adsorbing walks in three dimensions showing 
the four phases Desorbed-Extended (DE), Desorbed-Collapsed (DC), Adsorbed- Extended (AE) and 
Adsorbed-Collapsed (AC). Also shown are the dashed lines where simulations were performed. 



Pw = Pp = 0. Since Pp = the polymer is referred to as "Extended" and since Pw = the polymer 
is referred to as "Desorbed" . 

Now let us consider fixed = 1 and varying ujp. That is, we consider the collapse transition 
of a polymer fixed to a non- interacting wall as the intra-polymer attraction is increased. As the 
temperature is lowered (increasing the effect of the attraction) the collapse transition is understood 
to occur at the 0-point, where for high temperatures pp = (the Extended phase), while for low 
temperatures = 1/3 and Pp > (the Collapsed phase): this is expected to be a second order phase 
transition pLj. In both two and three dimensions this transition has been extensively studied (see 
for example [7] and references therein). 

On the other hand, if we consider fixed = 1 as the temperature is lowered {uJw is increased), 
the wall density is expected to change from pw = (Desorbed phase) for high temperatures to 
Pw > (Adsorbed phase) for low temperatures. This adsorption transition is also expected to be 
second order [2] and also has been well studied with this description being numerically well verified 
(again also see [7]). Note that for low temperatures the polymer is in a two-dimensional excluded 
volume state, where = = 3/4 [8|. We know that a two-dimensional polymer can also admit a 
collapse transition, so one naturally can conjecture a state where the polymer is both adsorbed with 
Pw > and two-dimensionally collapsed. A suitable two-dimensional order parameter to describe 
this transition is 

Up = hm ap{n) = lim „ . (5) 

n->oo n^oo Kg(n) 

When considering the full problem where both ujp and uJw can vary Vrbova and Whittington [Hill] 
conjectured, and numerically verified, a phase diagram with precisely these four phases: Desorbed- 
Extended (DE) phase with pp = and Pw = for small tOp and Uw, Desorbed-Collapsed (DC) 
phase with pp > and Pw = for large Up and (relatively) small w^;. Adsorbed- Extended (AE) 
with ap = and > for (relatively) small Up and large w^, and Adsorbed-Collapsed (AC) with 
cjp > and > for large ujp and ujw . 

Let us return to the bulk Desorbed-Collapsed phase where pp > and p^ = 0. One may now 
ask: how does p^ approach zero as the length is increased or, rather equivalently, how does the 
average number of wall visits (m^) scale as n is increased. It is believed (see below) that in the 
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Desorbed-Extended phase (m^) = 0(1), that is it is bounded. One may be tempted to assume that 
this is also the case in the Desorbed-Collapsed phase. However, recent work [9l [lOl [11] suggests 
that two different asymptotic behaviours can occur inside the latter. 

To understand this, let us consider the finite-sized extensive (reduced) free energy F{n;ujp,uJw) = 
— log(Zn). From thermodynamics we expect that 

F{n;iL!p,iL!^) = fh{iL!p,uj^u)n + o{n) , (6) 

and as stated above it is the analytic structure of the bulk free energy fb{ujp,ujw) that defines 
the phases of our model. Drawing on the standard polymer scaling in [12], the behaviour of the 
finite-size free energy in the DE phase is given by 

F{n; ivp, ujw) = hi^^p, iOw)n + (71 - 1) log(n) + 0(1) , (7) 

where 71 is an universal constant for the DE phase (an exponent, in fact). In the DC phase the 
scaling is expected to behave as 

F{n;u;p,ujyj) ~ /^(wp, u;^)n /s(wp, u;^)n^/^, (8) 

where we have a surface free energy fs{iOp,oOw)- This behaviour arises as the polymer assumes a 
dense liquid-like drop with a well defined surface that has area cx n^/^. We reiterate that this 
"surface" is not the wall but rather the surface of the liquid-like polymer drop. In the extended 
phase the polymer does not have a well defined surface and so = as in equation d?]). 
Importantly, the mean number of wall visits (m^„) can be calculated in the usual way as 

= jj:^!^^ . (9) 

Hence 

51og(cj^) d\og{ujy,) 

Since we assume that in both the Desorbed-Extended and Desorbed-Collapsed phases (m,^) = o(n), 
it is clear that fbi'^p^^w) should not be depend on cj^. Hence we have 

/ \ dfs{u}p,Ulw) 2/3 ^^^\ 

i'mw) TT7-7 — . (11) 

dlog(tj^) 

In the Desorbed-Extended phase fs{ujp, Uw) = 0, which derives (with some other weak assumptions) 
the {niiu) = 0(1) result for that phase. The question that then arises is the behaviour of fs{LOp,uJw) 
for the Desorbed-Collapsed phase. 

In [9l dni Hi] it was suggested that the Desorbed-Collapsed phase accommodates two different 
behaviours: one where /s(wp,u;^) depends on and so (m^) ~ n^/^, which was dubbed the 
"Surface Attached Globule" or SAG, and one in which fs{ojp,ujyj) is independent of w^, so that 
{rriyj) = 0(1) as in the Desorbed-Extended phase. We shall refer to this second situation as the 
"Fully Detached Globule". Physically the SAG phase can be pictured as the liquid drop partially 
wetting the wall, while the FDG is the state where the wall is dry. (The Adsorbed-Collapsed phase 
is the equivalent of a fully wet state.) This implies a surface phase transition with a singularity in 
fs{ujp,u!w) as iOw is varied, whereas fh{uip,u!w) is analytic at the same point and no bulk transition 
occurs. 

The evidence given in [911101 [TT] for the two 'surface' states was based on the analysis of relatively 
short exact enumeration data and some analysis of a directed walk model in two dimensions. It 
is therefore important to test this conjecture with data from longer walks. In this paper we have 
simulated, using a recently developed Monte Carlo algorithm [13], the model described above along 
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various lines in the parameter space for lengths up to 256. In a previous paper [14] we considered 
the bulk phase diagram and various low temperature finite size features. 

In [10] the boundary between the FDG and SAG was conjectured (see Figure [3] (a)). Using a 
zero temperature argument |T3], it is easy to see that for infinite cvp the transition from FDG 
to SAG occurs at uJw = 1- However, the question arises as to the form of the phase boundary 
for finite ujp. In [lOj it was conjectured that for finite ojp (but large enough so that the polymer 
is collapsed) the boundary occurs at values of w^^, greater than one. Hence, with small but finite 
attractive wall potential the polymer is expected to still be in the FDG state. On the other hand, 
if there is no entropic penalty to be paid for the globule to sit on the wall, this gives rise to an 
alternative hypothesis for the FDG-SAG boundary as shown in Figure [3] (b). We note that the 
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Figure 3: Hypothetical boundary for the surface phase transitions, (a) conjectured in |10| . (b) is 
an alternative that is suggested in this paper. 

phase boundary is unlikely to lie in the region for ujyj < 1; here the wall is repulsive and does not 
increase the surface entropy of the polymer, so one may argue that typical configurations have very 
few visits to the wall. Hence one conjectures that for uj^ < 1 the surface free energy fs{u}p,ujw) 
is independent of oj^ and (m^) = 0{1). The central question therefore is whether or not the 
polymer-drop partially wets the wall as soon as the wall potential becomes attractive. In the 
original hypothesis in Figured (a), a sufficiently attractive wall is required before the polymer-drop 
partially wets the wall. The alternative proposed in Figure [3] (b) conjectures that any attractive 
potential at the wall will induce the polymer-drop to partially wet the wall. 

So as to delineate the SAG-FDG phase boundary and to test for the existence of the SAG phase, 
we have studied three lines in the phase diagram (see Figure [2]) . One line at ujp = 2.0 was chosen 
since the 6 point is expected to be around LOp = 1.5 at the lengths considered here, and so by fixing 
iOp = 2.0 and varying uJw one explores the Desorbed-CoUapsed phase. We note that the position 
of the collapse transition should not move as is varied [1]. To further consider the difference 
between attractive and repulsive wall interactions on the Desorbed-Collapsed phase, we have also 
simulated along lines lOw = 0.8 (repulsive wall visits) and uJw = 1.3 (attractive wall visits). The 
second of these lines also allows us to search for the surface phase boundary proposed in Figure [3] (a). 
Note that the position of the adsorption transition for small ojp moves little as ujp is varied. 

In order to explore the Desorbed-Collapsed phase, we have used the FlatPERM Monte Carlo 
algorithm [13]. This algorithm estimates the density of states directly and so allows us to compute 
the partition function and related quantities for a wide range of LOp and lOs ■ The memory required to 
store the density of states grows with the cube of the length and so we were restricted to SAWS of 
maximum length 256. Each simulation sampled approximately 10^ confirmations. To reduce errors 
we performed 25 independent simulations and combined the resulting data. 
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Before we discuss the subtleties of the surface phenomena of the SAG and FDG phases we first 
verify that our simulations were correctly identifying the bulk phases. To do this we estimated the 
average end-to-end distance of the polymer which has been established to scale in the same manner 
as the (more computationally costly to estimate) radius of gyration. Hence we could estimate 
the size exponent v at various points. This confirms earlier work [1]. We considered the point 
{ujp,ijJw) = (1, 1), which is in the Desorbed-Extended phase, and so should have u 0.59: we find 
an estimate 0.60(2). Next, for the point {ujp^ujyo) = (1,2) which is expected to be in the Adsorbed- 
Extended phase and so should have v = 3/4 we find an estimate 0.75(1). Similarly, for the points 
{ujp,ijjw) = (2,0.8) and (2, 1.3) which are in the Desorbed-Collapsed phase (and we identify below 
to be in the FDG surface and SAG surface phases respectively) and so should have = 1/3 we find 
equal estimates 0.33(2). 

We begin our discussion of the surface phenomena by verifying the existence of the SAG phase. 
Let us consider the lines ujw = 1.3 (attractive wall visits) and oj^ = 0.8 (repulsive wall visits). 
In Figure m the mean number of visits m^t,(n; Wp, 1.3) divided by v?^^ is plotted against ujp for 
1.0 < ujp < 3.0 and n = 64, 91, 128, 181, and 256. The quantity m^{n) jr?!"^ would be expected to 
converge to a finite (non-zero) value only in the SAG phase. While there are clearly some corrections 
to scaling still evident at these lengths, rriw/m?/^ seems to be convergent to a non-zero value for 
ojp > 1.5, which is the rough location of the collapse transition at this length (the thermodynamic 
location has been estimated to be near 1.3). We have also checked that for ojp « 1.5 (that is, 
in the DE bulk phase) converges to a finite value. Hence the (finite-size location) bulk phase 
transition at 1.5 signifies a change in the length scaling of m^(n) from for small uip (DE bulk 
phase) to n^/^ for large ojp (DC bulk phase). 




1 1.5 2 2.5 3 



U}p 

Figure 4: The mean number of visits divided by v?^^ against ujp> 1 for fixed uJw = 1-3. 

In Figure[5l the raw mean number of visits raw[n] uJp, 0.8) is plotted against uJp for 1.0 < ujp < 3.0 
and n = 64, 91, 128, 181, and 256. As with uJu, = 1-3 one expects that for ujp < 1.5 (DE bulk phase) 
that raw{n) will converge to a finite value as n is increased and indeed this is the case. Additionally, 
the quantity myj{n) would be expected to converge to a finite (non-zero) value in the FDG phase. It 
appears that the collapse transition at around ojp = 1.5 only effects the value of limj^_>oo myj{n) but 
not whether it is finite or not. There is no indication of a change of length scaling when varying Wp. 

In Figure [H the variance of the number of visits is plotted against ujp for uj^ = 1-3 (a) and 

= 0.8 (b). In (a) a single peak in the variance is developing around Up = 1.4 as the length is 
increased, whereas in (b) no such peak is developing. We conclude a wall visit transition is occurring 
at around the same place the bulk contact transition occurs when lOw = 1-3 hut not when u)^ = 0.8. 

We now turn to considering analysis of the line LOp = 2.0. We first reinforce the analysis of 
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UJp 

Figure 5: The mean number of visits against ujp for fixed = 0.8. 
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Figure 6: Plot of the variance in the number of visits vs tUp for = 1.3 (a) and uj^j = 0.8 (b). This 
clearly shows a difference between SAG and FDG. 

Figure m and [5] by plotting the mean number of visits divided by along the line u>p = 2.0. We see 
that using 6 = we identify the FDG regime for uJw < 1 in Figure \7\ (a) . We also identify the SAG 
regime by using 5 = 2/3 in Figure [7| (b) for 1 < and at least tOw < 2, but probably for large lOw 
also. Finally the Adsorbed- Extended phase can be seen for > 3 where 5=1. 

Considering the variance in the number of contacts (Figure [8]), the bulk adsorption transition 
can be seen as a peak developing around uJw ~ 2.5, while another smaller peak can be identified just 
above 1.0 (see Figure [9]). While one might try to extrapolate the peak positions in Figure [9] to find 
the thermodynamic limit position, we find that our data is consistent with a wide range of limits 
ranging from 0.9 up to 1.15, depending on the choice of the correction exponent used to do the 
extrapolation. Certainly the peaks are moving to smaller values of tOy^ as the length is increased. 

We now return to the analysis of the line uJuj = 1.3. In Figure \T0\ we plot the variances in 
wall visits and bulk contacts for lengths 128 and 256 on the same figure, scaling the bulk variance 
so that the corresponding peak heights of the wall visit variances are roughly comparable. This 
enables us to make an easier comparison of the peak positions at the two lengths. We reiterate that 
only one peak occurs in the variance of the wall visits. If indeed this peak were associated with 
a surface transition occurring at a higher value of ujp than the bulk transition, one would expect 
to see the peaks of the variance of the visits extrapolating to a higher value of Up than the peaks 
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2 2.5 3 3.5 4 4.5 



(c) 

Figure 7: Plot of the mean number of visits divided by vs u>u! for ujp = 2.0. (a) shows the scahng 
in the FDG regime with 5 = 0. (b) shows the scaling in the SAG regime with 6 = 2/3. (c) shows 
the scaling in the thermodynamic adsorbed extended phase with 6 = 1. 
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Figure 8: The variance in the number of visits vs ujuu at ujp = 2.0 and length 256. One can clearly 
see two peaks indicating the possibility of two transitions. 
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Figure 9: The variance in the number of visits vs uJw at Up = 2.0 for different lengths around the 
SAG-FDG boundary. 



of the contacts variances. In both cases the peaks of the variances in wall visits occur at smaller 
values of ujp than the peaks of the bulk contacts, and both are moving to smaller iVp as the length 
is increased. We have investigated the scaling of the peak heights of the variances in wall visits and 
bulk contacts illustrated in Figure [TOl As expected the peak height of the (normalised) variance 
of the bulk contacts is consistent with a logarithmic scaling — an effective scaling as (log(n))^'^ is 
seen (one expects for very large n scaling of ((log(n))^/^^). The scaling of (normalised) variance 
of the surface visits is consistent with a logarithmic scaling using an effective form of (log(n))'''^ 
(or a weak power of -nP'^). Presumably, if both have logarithmic scaling this would be consistent 
with both peaks being associated with the bulk collapse transition. The behaviour observed does 
not support the scenario in Figure [3] (a) . This suggests that there exists one phase transition in 
the thermodynamic limit that gives rise to singularities in the scaling of the contacts and the wall 
visits. 




Figure 10: Variance in wall visits and bulk contacts at = 1.3 vs tOp at lengths 128 and 256. Note 
that the peak in the variance of the visits is to the left of the peak of the variance of the number 
of contacts. The variance of the number of contacts has been scaled so that the peak height of a 
given length is approximately the same as the peak height of the variance of the number of visits 
at the same length. 

Finally, we show typical examples of configurations of length n = 256 in Figure [11] at the points 
(upjUJii;) = (2,0.8) and (2,1.3) which we have argued are in the FDG and SAG surface phases 
respectively. One clearly sees a globule like conformation in both examples with a dense amorphous 
grouping of monomers. At (2, 0.8) only the monomer tethered to the wall lies in contact with the 
wall while at (2, 1.3) a fair proportion of the monomers of the surface of the globule lie on the wall 
as one would expect in a surface-attached globule. 
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(a) 



(b) 



Figure 11: Examples of typical configurations of length n = 256 at the points {ojp^ujw) = (2,0.8) 
(left) and (2,1.3) (right). 



We conclude by summarising that while we verify the existence of the SAG phase we find 
no indication of separate transitions other than occur in the bulk phase diagram or along zero 
interaction boundaries. This gives rise to a hypothesised surface phase diagram as illustrated in 
Figure [3] (b) . 
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